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Abstract. The injection rate of suprathermal protons into the diffusive shock acceleration process should vary 
strongly over the surface of supernova remnant shocks. These variations and the absolute value of the injection 
rate are investigated. In the simplest case, like for SN 1006, the shock can be approximated as being spherical in 
a uniform large-scale magnetic field. The injection rate depends strongly on the shock obliquity and diminishes as 
the angle between the ambient field and the shock normal increases. Therefore efficient particle injection, which 
leads to conversion of a significant fraction of the kinetic energy at a shock surface element, arises only in relatively 
small regions near the "poles", reducing the overall CR production. The sizes of these regions depend strongly on 
the random background field and the Alfven wave turbulence generated due to the CR streaming instability. For 
the cases of SN 1006 and Tycho's SNR they correspond to about 20, and for Cas A to between 10 and 20 percent 
Qf) , of the entire shock surface. In first approximation, the CR production rate, calculated under the assumption of 

spherical symmetry, has therefore to be renormalized by this factor, while the shock as such remains roughly 
spherical. 
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1. Introduction At the same time an essential physical factor which 

1 strongly influences the final CR acceleration efficiency, is 

; The time-dependent nonlinear kinetic theory of cosmic contained in our thcory as a free parameter. This is' the 

r> . ray (CR) acceleration in supernova remnants (SNRs) of ion injection rate; Le . the number of suprat hermal pro- 

CS Berezhko et al. CEIE|> and Berezhko & Volk l^!£0|2UUUJ, tons in j e cted into the acceleration process per unit area 

applied to the remnant of SN 1006^ Tycho's supernova and unit time It ig described by a dime nsionless injec- 

and Cas A (Berezhko et al. 2003 Volk et al. 12202 tion parame ter n that is a fixed small fraction of the 

Berezhko et al.|2003a|, has demonstrated that the exist- interstellar med i U m (ISM) particles entering the shock 

ing data are consistent with very efficient acceleration of fr(mt Assuming sphe rical symmetry, this injection implies 

CR nuclei at the SN shock wave, converting a significant a source term q _ Qj( r _ R ^ due to ( m0 no-energetic) 

fraction of the initial SNR energy content into CR energy. injection at the subshock m t h e diffusive transport equa- 

This energy is distributed between energetic protons and tion for the nudear CR distribution function f(r,p, t) 
electrons in a proportion similar to that of the Galactic 

— - = V(kV/) - w c Vf + — — p— + Q, (1) 



Recent Chandra observations (Bamba et al. [2303) dt 3 dp 

found the fine structure of the outer shock in SN 1006 
to be characterized by extremely small spatial scales of 



where Q s is written in the form 



the X-ray synchrotron emission. These structures agree Ni n jU\ 

very well with the above predictions and provide direct Q- s ~ 47rp? ^(p ~ Pinj)> (2) 



in] 



evidence for the efficient acceleration of nuclear CRs in 

this object (Berezhko et al. 2003b). with Ni n j = T)N\ and pi n j denoting the number of injected 

suprathermal particles from each unit volume intersecting 

Send offprint requests to: H.J.Volk the shock front and the momentum of the injected parti- 
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cles, respectively. Here r, t and p denote the radial coordi- 
nate, the time, and particle momentum, respectively; n is 
the CR diffusion coefficient; w c is the radial velocity of the 
scattering centers, u — V s — w is the flow velocity relative 
to the subshock at r = R s , V s = dR s /dt is the subshock 
velocity, N — p/m is the proton number density, and m 
denotes the particle (proton) mass. The subscripts 1(2) 
correspond to the upstream (downstream) region. 

According to nonlinear acceleration theory, the overall 
shock transition consists of a thin subshock whose thick- 
ness is of the order of the gyro radius of thermal ions 
heated in the shock compression, and a CR precursor 
whose much greater spatial extent corresponds to a mean 
diffusion scale of the accelerated particles. Ion injection is 
thought to occur at the subshock. When in the following 
we losely refer to "the shock" in the context of injection, 
then we always mean the subshock. 

Unfortunately there is no complete selfconsistent the- 
ory of a collisionless shock transition to the extent that 
it can predict the value of the injection rate and its de- 
pendence on the shock parameters for all directions of the 
shock normal relative to the external magnetic field vector. 
For the case of a purely parallel shock (where the shock 
normal is parallel to the external magnetic field) hybrid 
simulations predict quite a high ion injection rate (e.g. 
Scholer et al. lTMSl Bennett & Ellison HMty which corre- 
sponds to a value rj ~ 10 -2 of our injection parameter. 
Such a high injection is consistent with analytical theory 
(Malkov & VolkGSSlESa Malkov ITMSI) and confirmed 
by measurements near the Earth's bow shock (Trattner & 
Scholer HSU- 

We note however that in our spherically symmetric 
model these results can only be used with some important 
modification. In reality we deal with the evolution of the 
large-scale SN shock which expands into the ISM and its 
magnetic field. For example in the case of SN 1006 at the 
current evolutionary phase the shock has a size of several 
parsecs. On such a scale the unshocked interstellar mag- 
netic field can be considered as uniform since its random 
component is characterized by a much larger main scale 
of about 100 pc. Then our spherical shock is quasi-parallel 
in the polar regions and quasi-perpendicular in the equa- 
torial region. The leakage of suprathermal particles from 
the downstream region back upstream is to leading order 
dependent upon the shock obliqueness which is described 
by the angle between the ambient magnetic field direc- 
tion and the shock normal. It is most efficient for a purely 
parallel subshock and becomes progressively less efficient 
when the shock is more and more oblique (Ellison et al. 

Malkov & Volk nM)|> . Applied to a spherical shock 
in the uniform external magnetic field it would mean that 
only relatively small regions near the poles allow a suffi- 
ciently high ion injection rate which ultimately leads to 
the transformation of a significant part (more than a few 
percent) of the shock energy into CR energy, whereas the 
main part of the shock is an inefficient CR accelerator 
(Volk 20QTJ. In this case also the field amplification due 
to the CR streaming instability (Volk HnMI and in partic- 



ular Lucek & Bell 120001 Bell & Lucek l2UtJT)l occurs only 
near the poles, strongly amplifying the synchrotron emis- 
sivity there. Such a picture is consistent with the observed 
morphology of SN 1006: the most intense synchrotron ra- 
diation comes from two spots (Koyama et al. 1995; Allen 
et al. 2001), which we associate with the polar regions. 

This very simple picture probably holds only for 
"ideal" cases like SN 1006 which lies far above the Galactic 
plane in a very low density and apparently quite uniform 
environment. For other objects, presumably already for 
Tycho's SNR, and certainly for all core collapse SNRs, 
the situation is more complex. Yet the physical arguments 
which we shall use should apply to all of them with ap- 
propriate modifications. 

In this paper we quantitatively consider the system- 
atic variation of the ion injection rate across the SN shock 
surface in a simple approximation, taking the structure 
of the ambient magnetic field into account. For SN ex- 
plosions into a circumstellar medium that is not substan- 
tially modified by the mass loss from the progenitor star 
(i.e. SNe Type la, and core collapse SNe whose progenitor 
stars have zero age main sequence masses below about 15 
Mq), we shall demonstrate that the size of those shock 
regions, where efficient injection leading to efficient CR 
acceleration is expected to take place, depends strongly 
on the random background field and the Alfven wave tur- 
bulence generated by the CR streaming instability. For 
the case of SN 1006 this corresponds to about 20% of the 
entire shock surface, consistent with the observations. A 
similar result is obtained for Tycho's SNR. For the ex- 
treme case of Cas A, the final phases of stellar evolution - 
the red supergiant (RSG) and the subsequent Wolf-Rayet 
phase - have produced a complex circumstellar pattern of 
the magnetic field. It is characterized by an essentially az- 
imuthal mean field of stellar origin with superposed MHD 
waves, and strongly modified on large scales by instabili- 
ties due to the radial forces which have their origin in the 
late appearance of a high-luminosity Wolf-Rayet phase. 
An approximate calculation gives the result that only on 
a fraction of about 10 to 20 % of the shock surface ion 
acceleration proceeds efficiently. 

Electron injection is a rather different problem and 
our understanding is much poorer. In many ways, how- 
ever, electron injection can be discussed seperately. We 
will briefly adress it in section 6. While it is an essential 
ingredient for the synchrotron and the inverse Compton 
emission, to lowest order it neither affects the ion acceler- 
ation nor the overall SNR dynamics. 



2. Spherical shock in a uniform external magnetic 
field 

Fig. ri] schematically illustrates a spherical shock in a uni- 
form ambient magnetic field B. In the simplest MHD ap- 
proximation, the magnetic field structure in the down- 
stream region is determined by the compression of the 
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Fig. 1. Schematic form of the interstellar magnetic field 
lines, modified by a spherically expanding shock (thick cir- 
cle). 

component perpendicular to the shock normal and is de- 
scribed by the relations: 

^2\\ = Bi^, B 2 ± = (jBi±, (3) 

where a is the shock compression ratio, B\\ — BcosO, 
B ± =Bsm6. 

Since all the scales which characterize the motion of 
thermal particles are very much smaller than the shock 
size, a plane wave approach can be used for their de- 
scription. Fig. illustrates the magnetic field structure 
near the subshock front. The cold upstream plasma is ad- 
vected with speed u\ = V s towards the shock front, is com- 
pressed and heated, and flows with speed u 2 = u\ja into 
the downstream region. Very fast particles from the ther- 
malised downstream population, whose velocity exceeds 
some critical value Vi n j and which move towards the shock 
front, arc able to overtake it and to penetrate into the up- 
stream region. Since the particle mobility increases with 
their speed v, these particles can be considered as diffu- 
sive and they gain energy in their random walk across the 
shock. Asymptotically, for velocities v 3> Vi n j (e.g. Malkov 
& Volk 1995 ), their pitch angle distribution approximates 
an isotropic distribution. 

The number of these injected particles is determined 
by the structure of the subshock transition (Malkov|l998 
see Fig. 1 from Malkov & Volk 1998). The most important 
physical parameter which determines this number is the 
particle velocity parallel to the shock normal U|| : all par- 
ticles which have vu > Vi n j are considered as injected. It 
is assumed that the mean free path of the particles which 
have a speed higher than the threshold Vi n j exceeds the 
thickness of the subshock. 

Assuming the injected particles to come from the tail 
of a Maxwellian distribution we can write 

77|| = exp{-v? nj /v% 2 ), (4) 




Fig. 2. Schematic picture of the local magnetic field and 
the flow structure near the shock front. 

where vt 2 is the mean thermal speed of the downstream 
particle population. As mentioned above, numerical sim- 
ulations of parallel collisionless shocks give an expected 
injection rate m = r\(d\ = 0) ~ 10~ 2 (e.g. Scholer et al. 
HMR Bennett & Ellison HMty . that leads to the value 
Vi n j = 2vt2 of the injection velocity. 

We assume that the suprathermal particles are 
strongly magnetized. Therefore, for any magnetic field di- 
rection, the condition which selects injected particles is 
v \\ > Vi n j. Since v\\ = v cos 9 2 this means that only those 
particles are injected which move towards the shock front 
and have speed v > v- m j / cos 9 2 . Taking into account the 
relation 

cos 2 6» 2 = (1 + a 2 tan 2 6 1 y 1 (5) 
one finds 

v(ei)^vl +a2tan2ei - (6) 

According to this relation the injection rate goes down 
quickly with increasing upstream angle 9%, as illustrated 
by the curve which corresponds to SB = in Fig. [3J for 
Q\ = 10° the injection rate is one order of magnitude, and 
for Q\ = 20° it is four orders of magnitude smaller than 
for 6»i = 0°. 

This simple relation © is a direct consequence of our 
assumption that all particles are able to propagate only 
along the magnetic field lines. In reality they will to some 
extent also undergo cross field motion, either due to drift 
motions or due to cross field diffusion. On the other hand 
this cross field motion can be substantial only for a rel- 
atively high random field component. As shown below, a 
high-amplitude random field allows efficient particle injec- 
tion also within our simplified approach. 

It is important to note that there exists a so-called 
critical injection rate from the solution of the nonlinear 
kinetic Eq. coupled with the hydrodynamics of the 
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Fig. 3. The injection rate 77 as a function of the upstream 
angle 6\ between the ambient magnetic field and the shock 
normal, for different amplitudes SB/B of the upstream 
random field (solid lines), given that rj\\ = 10 -2 , for the 
current parameters of SN 1006. Dashed lines represent 
the critical injection rate n cr i tl the minimal injection rate 
Vmin required to provide significant random field ampli- 
fication, and the injection rate t]qcr which is equivalent 
to GCR reacceleration. The thick solid line represents the 
expected current injection rate; the vertical dotted line 
indicates the corresponding angular region of efficient in- 
jection/acceleration. 



thermal gas. The value of this critical injection rate is 
approximately determined by the expression (Berezhko & 
Ellison QUI! 



v \ - a / Pmax \ 
Vcrit = 10 

c V mc I 



-1/4 



(7) 



where c is the speed of light and p ma x is the maximum 
momentum of the accelerated CRs. It divides the region 
T) > Tjcrit of the efficient CR acceleration, when a signif- 
icant fraction of the shock energy goes into CR energy, 
from the region rj < rj cr it of inefficient CR acceleration. 
For V s = 3000 km/s and p m ax ~ 10 5 mc the critical injec- 
tion rate is r\ cr i t = 6x 10~ 5 . This means that in the case 
of SN 1006, if we do not take into account the hitherto 
disregarded random magnetic field component, efficient 
CR production is expected to occur only within two po- 
lar regions with 9\ < 9 max — 14°. This angular width is 
considerably smaller than the one observed (e.g. Koyama 
et al.HSSSl Allen et al.lMHT). 



3. Magnetic field fluctuations 

The existence of a random magnetic field component SB 
on scales large compared to the thickness of the subshock 



can change the value of the injection rate. To study this 
effect we assume that the ambient field 



B' 1 =B 1 + 5B 



(8) 



consists of two components, the uniform field B±, and 
a superimposed, isotropically distributed random compo- 
nent SB. If the spatial scale of the random component 
is much smaller than the shock size R s , one can find the 
mean injection rate by averaging over the directions of the 
random field: 



1 / t^-. l+tr" tarr" 

V=7Z disarm 



4tt 



'II 



(9) 



Here 9[ is the angle between B'^ and the shock normal 
direction n, in analogy to Fig. [3 The averaging proce- 
dure either corresponds to an average over a time interval 
which is short compared with the shock age and large 
compared with the relevant periods of resonantly scatter- 
ing field fluctuations for injected particles, or to a spatial 
average within the shock region whose size is small com- 
pared with the shock size and large compared with the 
scales of the relevant scattering field fluctuations. 

The averaged injection rate as a function of angle 6\ 
for different random field amplitudes SB/B is shown in 
Fig. El At small angles 9\ , where the shock is almost purely 
parallel, field fluctuations make it more oblique for some 
fraction of the time. The opposite is true for large 9\. 
Therefore, as one can see in Fig. |3 the existence of the 
random field component leads to a decrease of the injec- 
tion rate at small angles and to an increase at large angles 
01, so that at the highest values of the turbulent field 
5B/B — 1, compared to an ideal parallel shock, the in- 
jection rate is reduced by almost two orders of magnitude 
to a value 77(^1) ~ 10~ 4 which becomes almost uniform 
across an angular range 9\ < 63°. 

4. Selfconsistent turbulent field 

The random magnetic field component SB can be created 
selfconsistently by the CR streaming instability in the up- 
stream region rBell lT^75l Blandford & Ostriker rn37g|l . The 
expected amplitude of the Alfven waves excited due to the 
CR streaming instability is determined by the expression 
(McKenzie & Volk [EES Bell & Lucek lM)T]l 



5B_y = V s P c 
B 



c a pvr 

where c a is the Alfven speed, 
An 

T 



p. 



dppvf(p) 



(10) 



(11) 



is the CR pressure at the shock front, p in j and p max are 
the injection and maximum CR momenta, respectively, 
and / is the solution of Eq. (JTJ in which the injection rate 
appears as a given source term. Since we consider effects 
which can increase an initially very low injection rate, the 
shock modification by the CR pressure is neglected and a 
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plane wave approximation is used. Within this approach 
the CR distribution function at the shock front can be 
written in the form (e.g. Berezhko et al. 1996) 



j, _ qNinj f P V 



(12) 



where q = 3ui/ (u± — u%). Taking into account that in the 
case of a strong unmodified shock (q = 4) the appropriate 
value of the injected particle speed is v in j ss 2V S , and that 
relativistic particles with mc < p < p m ax provide the main 
contribution to the CR pressure, we can write 



SB \ _ 8aq p max 
B I 3c a mc 



(13) 



The quantity hi(p max /mc) depends at most logarithmi- 
cally on (SB/B) 2 . As a consequence it can be seen from 
this relation that, once the selfconsistent Alfven wave field 
SB exceeds the background ISM fluctuation field 8Bo (see 
Sect. 7), then any initially low injection rate leads to the 
growth of the random magnetic field in the upstream re- 
gion which in turn leads to an increase of the injection 
rate (see Fig. Equating (SB/B) 2 to the background 
level (SB/B) 2 , we can find the minimal initial injection 
rate 

3c a (SB\ 2 



8c\n(p max / mc) \B J Q 

On account of the Alfven wave generation by the CR 
streaming instability the expected high efficiency injec- 
tion region is bounded by the polar angle 9 ma x determined 
from the relation 



v (SB,t 



(15) 



where rj(SB, 9\) is the function shown in Fig.|3for a given 
random field amplitude 5B. Within the region 9\ < 9 max 
the initial injection rate is high enough so that acceler- 
ated particles are able to increase the level of the tur- 
bulence level which in turn increases the injection rate. 
This selfconsistent nonlinear amplification also increases 
the mean magnetic field strength to an effective mean field 
B whose difference to the actual field defines an effective 
SB (Lucek & Bell 2000 ). The process can be assumed to 
end when the amplitudes of the Alfven waves become so 
high that SB ~ B for the effective quantities and their 
further growth is prevented by strong nonlinear dissipa- 
tion processes. Identifying therefore (SB/B) in Fig. [31 for 
9\ < 9 max with the ratio of the effective quantities, the 
expected selfconsistent injection rate corresponds to the 
curve T](SB = B,9i). 

The amplification of the field to an effective field large 
compared to the external field Bi is an important as- 
pect of CR production, since it determines the value of 
the maximum CR energy. At the same time, this strongly 
modified upstream field remains completely randomized. 
Therefore we assert that the injection rate within the re- 
gion 9\ < 9 max is not sensitive to the specific value of the 
upstream magnetic field B±. 



5. Background cosmic ray acceleration 

Suprathermal particle leakage from downstream is not the 
only mechanism that supplies particles to the diffusive 
shock acceleration process. Galactic cosmic rays (GCRs) 
are also subject to further acceleration, sometimes also 
called re-acceleration. We assume that the majority of 
GCRs which are participating in this re-acceleration have 
a momentum p ~ mc and therefore their number density 
can be estimated as 



Ngcr = eGCR/(mc 2 ), 



(16) 



where cqcr. ~ 0.6 eV/cm 3 is the GCR energy density. 
Comparing the number of CRs with momenta p > mc, 
produced in shock acceleration at a given injection rate 
m with Ngcr-, we can define a minimal injection rate of 
suprathermal protons 



Vgcr 



cNgcr 
8V S N H ' 



(17) 



below which GCR re-acceleration becomes more efficient 
than acceleration of suprathermal particles. 

It can be seen from the above expression that the role 
of GCR re-acceleration increases with decreasing shock 
speed and/or ISM gas number density. Therefore it can 
become significant in the diluted ISM during late SNR 
evolutionary phases, or for SN explosions in the hot and 
low density ISM of elliptical galaxies (Dorfi & Volk lTM)}) . 

Since re-acceleration of GCRs is almost independent of 
the polar angle 9\ fDrurv H983|) . efficient CR production 
will occur practically over the whole SN shock surface, if 
rjGCR. exceeds the critical injection rate 

In the case t\gcr < Vera the role of GCRs can be still 
significant if i]gcr. > Vmin- Since at large angles 9\ GCR 
reacceleration dominates over the suprathermal particle 
acceleration, Alfven wave generation due to reaccelerated 
GCRs is also higher: they produce in the upstream region 
an Alfven wave field 



(18) 



Therefore, efficient injection in this case occurs within the 
angular range 9\ < 9 max where 9 max is determined from 
the relation rj(SB, 9 max ) = VGCR, with SB from expression 

6. Electron injection 

The injection of electrons is a different problem and much 
less well understood physically. We shall give here a very 
brief discussion in order to connect synchrotron and, pos- 
sibly, inverse Compton emission to the dynamics and ra- 
diation properties of accelerated nuclei in SNRs. 

Suprathermal electrons from the hot downstream re- 
gion cannot resonantly scatter on the MHD waves pro- 
duced by the escaping ions, and other wave types have 
to be investigated (Levinson 1996). On the other hand, 
in more or less perpendicular shocks reflection of part 



(s: 


\ 2 _ VGCR 






Vmin 
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of the incoming ion population back into the upstream 
plasma is possible, exciting electrostatic waves there in 
which electrons can be energized. This energization can 
certainly inject these electrons into the diffusive shock ac- 
celeration process if the electrons reach GeV energies from 
which point on they accelerate like relativistic ions of com- 
parable kinetic energies. Models of this kind have been 
investigated after the pioneering study of Galeev (1984) 
by Galeev et al. (1995), McClements et al. (1997) and 
Dieckmann et al. (2000). This suggests that, in contrast 
to ions, electrons may be even best injected at rather per- 
pendicular shocks. Since the spatial and temporal scales of 
this initial electron energization are very short, this argu- 
ment can presumably be applied locally to any part of the 
shock surface where the shock normal is highly oblique 
to the local instantaneous magnetic field. Since further- 
more in the quasi-parallel shock regions the energetically 
dominant accelerating ions create large amplitude MHD 
waves, there may thus be parasitic electron injection and 
acceleration also in these regions. Altogether, electron in- 
jection might therefore occur more or less everywhere over 
the SNR shock surface. Whether this can hold also for 
subsequent acceleration to the highest energies, is a quite 
different matter, and less than clear. 

Let us nevertheless assume that electrons are uni- 
formly injected and even accelerated across the SNR 
shock. The resulting synchrotron emission is oc 
and will therefore still occur predominantly in the polar 
regions with their strong turbulent field amplification, de- 
spite the average field compression in the equatorial re- 
gion (see next section). The Inverse Compton gamma-ray 
emission may, however, be fairly uniform over the SNR 
surface. 

7. Results and discussion 

As a representative case, we first apply the above for- 
malism to SN 1006. For the relevant SN 1006 param- 
eters V s = 3200 km/s, and Nh — 0.1 cm" 3 we have 

VGCR = 7 X 10" 8 . 

Since the interstellar random magnetic field is dis- 
tributed over a wide range of scales A, only part of this 
spectrum with scales A < X rnax has to be considered as 
the small scale field. The upper scale can be taken as the 
diffusive length of the highest energy CRs accelerated, 
i.e. X m ax = l(pmax)- According to our numerical results 
(Berezhko et al. 2002), for the current evolutionary phase 
of SN 1006, p max = 4 x 10 5 mc and l(p max ) = 0.08 R s . 
Assuming that the energy density per unit logarithmic 
scale interval of the background ISM turbulent field as 
function of the spatial length A increases according to the 
Kolmogorov law 

< i9 > 

we find that for scales smaller than X max = 0.08 R s it 
has a value [5B/B)l = f* max dlnXE w /(B%/8w) = 5.3 x 



10 -2 , or (5B/B)o — 0.23, taking into account that the 
main scale is L = 100 pc. According to Eq. (fbEjl . for a 
typical ISM magnetic field B = 3 fiG and a moderate 
maximum CR momentum p max — 10 3 mc, we have in this 
case rimin = 1.7 x 10~ 7 . Since r) GCR < r) min , GCR re- 
acceleration does not play an important role in this case. 

As one can see from Fig. 01 the line f](9i) which cor- 
responds to SB/B — 0.23 intersects the minimal level 
Vmin = 1-7 x 10~ 7 at 8 max — 31°. This means that an 
initial injection level represented by this curve rj{9\) leads 
within the whole region 9\ < 31° to progressive growth 
of Alfven waves, which in turn leads to a corresponding 
increase of the injection rate. As argued above, this posi- 
tive backreaction drives the field amplification to the point 
SB/B ~ 1. Therefore the expected injection rate r\ « 10~ 4 
at 6*i < 31° corresponds to the curve rj(SB = B,6i). Since 
f] > Vcrit, efficient CR production will occur in this angu- 
lar region. For 9\ > 9 max the expected injection rate 77(6*1) 
goes along the curve rj(6B/B = 0.23, #1) and then along 
V = VGCR- For all 9i > 9 max the expected injection rate 
is lower than the critical rate. 

Then we have the situation where a sharp boundary 
9i = 9 max separates the regions of efficient and inefficient 
CR injection/acceleration. One has to expect that in re- 
ality this boundary is smoothed as a result of some addi- 
tional physical process. Possibly the most important fac- 
tor is the cross-field diffusion of CRs. Due to the high level 
of selfconsistent turbulence within the region 9\ < 9 max , 
CR diffusion is almost isotropic. High energy CRs with 
p > mc will therefore be able to penetrate diffusively 
through the boundary 9i — 9 max in the upstream region. 
They can also be accelerated in a finite region 9\ > 9 max . 
An approximate CR diffusion length across the regular 
magnetic field in the upstream region is their parallel 
diffusion length l(p max ). It corresponds to the angle in- 
terval A6>i = (l/R s ) rad w 5°. This means that the 
smoothed region of efficient CR acceleration extends up 
t° Cm = 9 max + A9 « 36°. Therefore efficient particle 
injection/acceleration is expected to occur within a bipo- 
lar region of about 20% of the shock surface. Its size cor- 
responds rather well to that of the observed bright X-ray 
synchrotron emission regions of SN 1006 (e.g. Allen et al. 
I20UTJ) 1 . 

7.1. Renormalization 

According to the above estimate a substantial part of the 
shock still efficiently injects and accelerates CRs. In ad- 
dition, the overall conservation equations ensure an ap- 
proximately spherical character of the overall dynamics. 

1 We shall not go into speculations here whether the non- 
linear wave amplification process can spread into SNR surface 
regions where the shock is practically perpendicular, so to say 
by itself. While such a scenario can not be excluded theoret- 
ically at this time, we believe that it is more fruitful to ask 
whether the observations give any hint in this direction. From 
present knowledge they do not. 
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Therefore, we assume the spherically symmetric approach 
for the nonlinear particle acceleration process to be ap- 
proximately valid in those shock regions where injection 
is efficient. To take this injection fraction 

/ re = l-COsCa* (20) 

into account, we need then to introduce a renormaliza- 
tion factor for the overall nuclear CR acceleration effi- 
ciency, and for all the effects which it produces in the 
SNR. According to the above estimate its value in the 
case of SN 1006 is / re 0.2. 

7.2. Synchrotron emission 

The amplification of the magnetic field in the region of ef- 
ficient acceleration leads to an effective downstream field 
strength of 120/xG strength that exceeds the mean inter- 
stellar field of about 3/iG by a factor of 40, according to the 
comparison of the synchrotron measurements for SN 1006 
with the theoretical model (Berezhko et al. 12002(1 . Outside 
the ion injection region, in the "equatorial" region, there 
is by definition no shock modification and the interstellar 
field is at most compressed by a factor of 4, correspond- 
ing to a strong adiabatic shock in the thermal plasma. 
In fact, the field compression factor is on average smaller 
than 4 because the shock is truly perpendicular only on 
the equator itself. Therefore the enhancement of the local 
synchrotron emissivity in the polar regions relative to the 
equatorial region exceeds 10 15 w 32 for the hard parti- 
cle spectrum oc p" 2 . The spatially integrated emissivities 
must be adjusted by the ratio / re /(l — /re) = 1/4 of the 
polar to the equatorial areas. Therefore the global excess 
emissivity of the polar regions is expected to be 8 times 
larger than the equatorial emissivity on account of the 
higher magnetic field. As a consequence we expect the 
polar regions to dominate also the synchrotron radiation, 
in contrast to the arguments of Ratkiewicz et al.(1994) 
who assumed that the uniform ISM field is merely MHD- 
compressed by the spherical shock, which would then obvi- 
ously lead to maximum emissivity in the equatorial region. 

7.3. The case of Tycho's Supernova 

Until now we have given numbers which refer to SN 1006. 
Also the discussion of morphological effects was done for 
this ideal case. Even though Tycho's SN was also of Type 
la, its present X-ray emission is dominated by line radi- 
ation, and not by nonthermal emission, like SN 1006. In 
addition, the ambient interstellar medium around Tycho 
appears to be unexpectedly nonuniform, cf. Reynoso et al. 
(1997). As a consequence, it is not clear whether we can 
expect a nonthermal X-ray morphology as simple as that 
of SN 1006. Nevertheless, the magnetic field topology is 
analogous and therefore we shall calculate the renormal- 
ization factor in an analogous form. 

With a present radius R s rs 2.7 pc and a mean ambient 
density N H = 0.5 cm" 3 , Tycho is smaller than SN 1006 



and the ambient density is considerably higher, whereas 
the present shock velocity V s = 3100 km/s is about the 
same. This gives r\ m %n = 2.4 x 10~ 6 (SB / B)q. Similarly, we 
have {8B/B)l = 3/2(A maa; /Lo) 2/3 with X max = 0.27 pc. 
Taking the turbulent main scale Lq proportional to the 
critical SNR radius at cooling which in turn is propor- 

— 1/3 

tional to N H , we have Lq = 58.5 pc for Tycho, if 
L Q = 100 pc for SN 1006. This gives (SB/B) = 0.2, 
and finally f re = 0.18 for Tycho, roughly similar to the 
value for SN 1006, as used in Volk et al. <|2UU2)) . 

7.4. Late phases 

Since the energy density of the random component of the 
background field (5B/B)q oc R 2 J 3 increases during the 
SNR evolution, the size of the region of efficient injec- 
tion becomes progressively larger as well. In the extreme 
case of a shock with size R s ~ 100 pc - provided it is 
still strong enough - efficient injection/acceleration oc- 
curs across the entire shock surface due to the completely 
randomized background field on the scale R s . This situa- 
tion makes the contribution of the late SNR evolutionary 
phases more important, and therefore the resultant CR 
energy spectrum produced in such SNRs is expected to 
be steeper compared with the model prediction in spheri- 
cal symmetry. 

According to expression (JZJ) the critical injection rate 
goes down proportionally to the shock speed V s , whereas 
the injection rate (|17|l due to GCRs is inversely propor- 
tional to V s . Therefore the significance of GCR reacceler- 
ation in the sense of an injection mechanism progressively 
increases during SNR evolution. When the shock speed 
drops to the value V s rs 35/-\/A#/(l cm" 3 ) km/s, tjgcr 
exceeds r\ cr it and then the GCRs alone lead to efficient 
CR acceleration across the whole shock surface. Since the 
GCR chemical composition differs from the ISM composi- 
tion, GCR re-acceleration can play a role in determining 
the resultant chemical composition. 

7.5. Renormalization factor for Cas A 

A much more complicated situation arises in SNRs where 
the circumstellar medium is strongly influenced by the 
wind from the progenitor star. Cas A is a prominent exam- 
ple. According to the analysis of the thermal X-ray emis- 
sion, and consistent with the observed overall dynamics 
of Cas A, the supernova shock expands into an inhomoge- 
neous circumstellar medium strongly modified by the in- 
tense wind of the progenitor star (Borkowski et al. ,1996 ) . 
It consists of a tenuous inner bubble, created by the Wolf- 
Rayet (WR) wind, a dense shell of swept-up, slow red 
supergiant wind (RSG) material, and a subsequent free 
RSG wind. The ambient magnetic field structure is not 
well known in this case. Therefore it is not possible to per- 
form an equally definite analysis of the expected injection 
rate as for Type la SNe. As we shall demonstrate below, 
quite a reasonable estimate can nevertheless be obtained. 
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Fig. 4. The renormalization factor / re as a function of the 
random large-scale magnetic field component, for different 
values of the small-scale field component SB S . 

Due to the progenitor star's rotation, the mean large- 
scale magnetic field is expected to be almost purely tan- 
gential to the shock surface in the wind material. At the 
same time, as shown in numerical simulations (Garcia- 
Segura et al. 199^, the interaction of the fast WR with 
the slow, massive RSG wind and the shell formation is 
accompanied by a long-wavelength instability. We assume 
that this instability is able to strongly randomize the pre- 
existing field in the shell, and to possibly even amplify it. 
Therefore the ambient circumstellar magnetic field seen by 
the SNR shock in the shell should consist of three compo- 
nents 

B 1 = (B)+5Bi+SB s , (21) 

a mean field (B) which is practically everywhere tangen- 
tial to the shock surface, a large-scale random component 
(SBi), and a small-scale random component (SB S ). The 
small-scale field component consists of the copious amount 
of MHD waves, essentially Alfven waves, emitted by the 
central RSG, and is expected to follow a nonlinear cascade 
towards short scales, taken here to be of the Kolmogorov 
type. 

We shall assume an analogous magnetic field struc- 
ture in the outer free RSG wind region (Berezhko et al. 
20033). The relevant instability there is that of radiation 
pressure-driven winds (Lucy &; White 119801 for a more 
recent discussion, see Lucv ll98"4l which we consider here 
in terms of the dynamical effect of the intense WR star 
radiation field on the free RSG wind. 

During the SNR shock propagation through the RSG 
wind material, the variation of the random magnetic field 
component will lead in both phases - compressed shell 
and free RSG wind - to a number of irregular regions 
(spots) on the shock surface where the shock is quasi- 



parallel, whereas the main fraction of the shock is quasi- 
perpendicular 2 . 

The fraction of the shock surface f re where the shock 
is quasi-parallel and where therefore efficient injection oc- 
curs, depends on the ratios of the values of the magnetic 
field components (B), 5Bi and 8B S . In order to estimate 
the expected value of the renormalization factor f re for 
given values of (B), 5Bi and SB S , we first calculate the 
injection rate 77(^1 ) by performing the averaging proce- 
dure over all possible directions of the small scale random 
field B s according to the expression where in this 
case B± = (B) + Bi and SB = SB S . Subsequently, we 
can calculate the value of the solid angle Af_ which in- 
cludes those components of the large scale random field 
8B1 whose directions lead to efficient injection, that is 
77 (#1) > r\ m in ■ In the spirit of Fig. 3, we then assume that 
above r\ m in injection rises nonlinearly as a result of strong 
acceleration, if the typical scales of the high-injection flux 
tubes exceed the size l(p ma x) ~ 0.1i? s of the CR precur- 
sor. This is expected to be the case at least in the com- 
pressed shell of Cas A, judging from the numerical results 
of Garcia-Segura et al.(1996). In the outer free RSG re- 
gion we have also assumed this to be true in the model of 
Berezhko et al. I|2003afl . However, it is not clear how large 
the radial scales of the radiation-induced instability are, 
even though the scales of the frayed shock due to the fast 
moving ejecta knots should be the same as in the shell. 
Therefore the efficient acceleration also in the free RSG 
wind remains an assumption which we shall also make 
here. Quantitatively, the radio and X-ray flux from the 
dense shell dominates the overall nonthermal emission in 
these wavelength regions at the present epoch (Berezhko 
et al. 12003a). so that the contribution from the free RSG 
wind is not a decisive factor in the interpretation of the 
gamma-ray emission from Cas A. 

Since the gas number density in the shell is N g = 
10 cm~ 3 and the magnetic field value is about Bq = 
200 fj,G, the Alfven speed is c a = 130 km/s which gives 
r] min = 1.4 x 10- 5 (SB/B)l, where SB = SB S and B = 
(B) + SBi. Because the large-scale field component 6B1 
is randomly distributed, the value of the renormalization 
factor is determined by the simple expression 

f re = Afi/(47r). (22) 

The value of the renormalization factor f re as a function 
of 5Bi/(B), calculated for different ratios of small-scale to 
large-scale field components SB s /5Bi, is shown in Fig.QJ 
One can see that for SBi ~ (B) the renormalization fac- 
tor f re has values between 0.1 and 0.2 if the small-scale 
random component has a relative amplitude SB s /SBi be- 
tween 0.2 and 0.4. Such relative amplitudes appear quite 
realistic. 

2 Similar regions with irregular shock normal angles relative 
to the magnetic field will be produced by the clumpy nature 
of the SN ejecta - the fast moving knots - which give the SNR 
shock a rather frayed appearance. For lack of better knowl- 
edge we shall lump the two effects together here into a single 
component 8B1. 
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We conclude that the value f re = 0.15, which yields 
a satisfactory description of all relevant properties of the 
emission generated in Cas A by accelerated CRs (Berezhko 
et al. l2003a|) . is indeed consistent with the expected struc- 
tural properties of the circumstellar magnetic field. 

8. Summary 

The injection rate of suprathermal particles into the shock 
acceleration process depends strongly on the shock obliq- 
uity and diminishes as the angle between the ambient field 
and the shock normal increases. For the ideal case of a SN 
explosion into a uniform interstellar medium with a uni- 
form magnetic field B, efficient particle injection, leading 
to the conversion of a significant fraction of the kinetic en- 
ergy at a shock surface element, only occurs in relatively 
small regions near the "poles", reducing the overall CR 
production. The sizes of these regions depend strongly on 
the random background field and on the Alfven wave tur- 
bulence generated upstream of the shock due to the CR 
streaming instability. 

For SN 1006 which appears to approximate this ideal 
case, efficient CR production is expected to arise within 
two polar regions, where the SN shock is quasi-parallel. 
The relative size of these regions depends decisively on 
the amplitude of the random background field component 
SB: it changes from a fraction of about 0.07 to one of 
0.7 of the shock surface when the interstellar random field 
amplitude varies from SB/B = to 1, respectively. It is 
argued that the nonlinear backreaction of wave production 
on the injection rate leads to a definite size of the injection 
fraction. 

It is also argued that the actual spatially integrated 
nuclear gamma-ray emission from such objects can be ob- 
tained through a renormalization of the spherically sym- 
metric result by the same factor. 

For SN 1006, the calculated size of the efficient CR pro- 
duction regions which amounts to about 20% of the shock 
surface corresponds very well to the observed sizes of the 
bright X-ray synchrotron emission regions (e.g. Koyama et 
al.©95 Allen et al. l2001J) . This implies a renormalization 
factor of 0.2. 

In the case of Tycho's SNR, with its unexpectedly 
nonuniform ambient interstellar medium, we can not ex- 
pect such a simple X-ray morphology. However, at least 
the topology of the field should be the same. Therefore 
the renormalization factor f re has been calculated in an 
analogous manner. Although the relevant parameters dif- 
fer from those of SN 1006, we obtain the similar value 

fre = 0.18. 

The total size of the regions where efficient injection of 
suprathermal particles occurs and the relevance of GCR 
reacceleration are expected to increase during SNR evo- 
lution. This leads to a steepening of the resultant energy 
spectrum of CRs produced in SNRs compared with the 
spherically symmetric case. 

In the case of a circumstellar medium which is strongly 
perturbed by the mass loss of the progenitor star, like it 



exists for Cas A, efficient injection presumably takes place 
in a number of randomly distributed portions of the shock 
surface with quasi-parallel magnetic field. The fraction of 
the shock surface covered by these spots depends on the 
relative strength of the mean magnetic field, which is as- 
sumed to be tangential to the shock surface, and upon 
the amplitude ratio of the small-scale and the large-scale 
scale random field components. It is larger for lower mean 
field values and for higher small-scale random fields. It was 
demonstrated that quite reasonable magnetic field param- 
eters are consistent with efficient injection on about 15% 
of the shock surface, the fraction which is required in order 
to reproduce the observed properties of Cas A (Berezhko 
et al. l2008aj) . 
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